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1. Introduction 



Semiclassical approximation is a powerful and perhaps the most commonly used ap- 
proach to quantum mechanics for exploiting classical mechanics based on the idea that 
the former may be realized by supplementing the latter properly (see, e.g., [1]). In the 
path-integral language, it asserts that the transition amplitude between two arbitrarily 
given points may be approximated by summing up fluctuations around the classical path 
connecting the two points. Here, the existence of the classical path is assumed, not guar- 
anteed, in the first place. 

Caustics occur when this assumption breaks down. More precisely, when a family of 
classical trajectories focuses, the envelope of the trajectories forms a focal region called 
caustics. In one dimension, the region becomes a focal point, and hence a classical path 
connecting generic initial and final points a, b G IR exists if and only if b is the focal point 
specified by a. Given the action / of the system, this happens when the second variation 
of the action 5 2 I along a classical path starting from a vanishes, and as such their analysis 
is purely classical and constitutes a branch of the Sturm-Liouville problem. In particular, 
if we confine ourselves to the action I[x] = J dt L for a finite time interval [0, T] with the 
quadratic Lagrangian 3 

L=-x 2 -^X(t)x 2 -fi(t)x, (1.1) 

then 5 2 I = is equivalent to the condition that a solution u obeying the homogeneous 
equation of motion u{t) + X(t) u{t) = with u(0) = vanish at t = T. 

In quantum mechanics, the semiclassical approximation for the quadratic system is 
known to be exact, and in terms of the Morse index m(A) and the action I[x c \] evaluated 
for a classical path x c \(t), the amplitude for the transition between the two points a, b 
reads [2, 3] 

(L2) 

Thus the amplitude becomes singular when the harmonic potential A admits the solution 
u to become u(T) = 0, that is, when caustics occur. The appearance of the singularity 
suggests that, on caustics, a generalized prescription is required for semiclassical approxi- 
mation to treat the cases where the classical path does not exist. Such a prescription has 
been devised in the path-integral framework [4, 5] (cf. section 3.1 of this paper). Strictly 

Note that any Lagrangian at most quadratic in position and velocity can be brought into this form 
by partial integration. In this paper we use the dot to denote time derivative x = dx/dt, and put both the 
mass m of the particle and the Planck constant h (except section 3.2) unity for convenience. 
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speaking, the singularity does not arise in actual physical processes, because there usually 
exist higher order terms (cubic, quartic etc. in x) in the Lagrangian, or even if the quadratic 
form (1.1) provides a reasonably accurate description of the system, the potential A may 
not allow the solution u to vanish at t = T, which is generically the case. Nonetheless, 
physical phenomena pertaining to caustics do arise, and the quantum analysis for such 
phenomena occurring in the presence of higher order terms has been carried out in the 
path-integral by Schulman [4] and DeWitt-Morette [6] (see also [5, 7, 8, 9]; for physical 
aspects of caustic phenomena in wave theory and the relation to Thorn's theorem, see, 
e.g., [10]). 

In this paper, we wish to provide a study of quantum caustics characteristic to the 
quadratic Lagrangian (1.1) for a generic potential A and an external force \x. The aim of 
the study is three-fold: (i) to present a basic but self-contained result on classical caustics 
which can be used for a fuller analysis of caustics for general Lagrangians (because caustics 
are characterized by the quadratic part of the Lagrangian) , (ii) to obtain a formula for the 
transition kernel K(b, T; a, 0) which covers the case of caustics and introduce the notion 
of a quantum Jacobi field to discuss the physical implications of caustics at the quantum 
level, and (iii) to apply the result to examine Jevicki's correspondence hypothesis, which 
we now explain below. 

In an attempt to resolve certain discrepancies between instantons in QCD and the 
1/N expansion, it was argued by Witten [11] in the two dimensional <LP N nonlinear sigma 
model (which is a prototype of QCD) that classical instantons are eliminated by quantum 
fluctuations and disappear at the quantum level. To this assertion Jevicki [12] contended 
- under a few but crucial assumptions — that instantons do not disappear but show up 
in the form of simple poles rather than stationary points at the quantum level. His point 
was illustrated by the example of the quantum mechanical rotor, i.e., a free particle on 
the circle S 1 , where the correspondence 

stationary points •<=>- simple poles (1.3) 

may become transparent. Embedded in the plane IR , the rotor can be transformed into a 
two dimensional harmonic oscillator with an arbitrary potential A, where our result for the 
quadratic Lagrangian is applicable. Our detailed analysis will show that, even in this toy 
model, the assumptions made in [12] cannot be justified and, therefore, the correspondence 
hypothesis (1.3) cannot be sustained. In short, there is no firm reason for rejecting Witten's 
assertion. 
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The plan of the paper is as follows. In section 2 we provide a full account of classical 
caustics associated with the quadratic Lagrangian (1.1). To render the classical results 
most convenient for later use, we present them in a mathematical style. Readers who 
wish to see only the discussion of quantum caustics may skip this section. Section 3 is 
devoted to quantum caustics, where we first derive the kernel formula which admits the 
case caustics and then discuss the quantum effect associated with caustics by a simple 
Gaussian slit (gedanken-) experiment. Jevicki's correspondence hypothesis is examined in 
section 4. Section 5 contains our conclusion and discussions. An appendix is provided at 
the end to supplement the argument of section 3. 
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2. Classical Caustics 



In this section we review and summarize various aspects of the classical motion of 
a point particle governed by a quadratic Lagrangian. Our discussions focus especially 
on the occurrence of caustics and its consequences in the framework of classical physics. 
Part of the material presented is well-known, see e.g., [5], but is included here together 
with numerous other facets of the issue. In this way the following discussions illustrate 
the phenomenon of caustics in classical physics as opposed to that in quantum realm in 
later sections and further furnish a mathematical basis of the discussions in the following 
sections on quantum mechanics on S 1 . First we consider the Jacobi equation, which in 
the case of quadratic Lagrangian is the homogeneous part of the equation of motion, and 
discuss some general properties of its solutions. Next we concentrate on the phenomenon 
of caustics and discuss two of its main characteristics, the one being the constant stretching 
factor between the final and the initial points, and the other being the Morse index. 

2.1. Quadratic Lagrangian and Jacobi equation 

Let us consider a point particle governed by the quadratic Lagrangian (1.1), whose 
equation of motion reads 

d 2 



A\x(t) = /j(t) , where A\ := — 



(2.1) 



For the discussion of caustics it is useful to consider the Jacobi field J(p,t) (see e.g., [5]) 
which describes the spread of paths from an initial point when varying the momentum. 
Given a one-parameter family of classical paths x(p, t) characterized by the initial momenta 
x{p, t) = p at t = to, it is defined by 

J M (2.2) 

It satisfies J(p, to) = by definition and solves the Jacobi equation, which in case of the 
quadratic Lagrangian (1.1) is the homogeneous part of (2.1). Hence, the Jacobi equation 
can itself be interpreted as the force-free equation of motion, 

A x u(t)=0, (2.3) 

and as such contains a wider class of solutions (not necessarily zero at t = 0) which may 
be added to a solution curve of the full equation of motion (2.1) to yield other extremum 
paths. The investigation of the solution curves to (2.3) not only leads to a characterization 
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of caustics via the Jacobi field itself, but reveals further properties of the caustics, see 
in this context also [6]. In particular, the relationships between two linearly independent 
solutions of the homogeneous differential equation of second order (2.3) help us to clarify 
the notion of the stretching factor and the index which will be introduced later on. 

As is well known, given initial data x(to) = xq and x(to) = xq at to, the differential 
equation (2.3) possesses a unique solution. If A(t) is continuous, the solution is of class C 2 , 
i.e., has continuous derivatives of order two, and when A(t) is only piecewise continuous, 
the solution curve is at least C 1 . In what follows, we shall use frequently the following 

Lemma 1. Let v(t) and u(t) be solutions of the differential equation (2.3), and let t\ and 
ti be two times. Then we have 



vu 



= vu 

tl 



(2.4) 



For the proof, we observe from (2.3) that 



vu 



** f 2 . d , , r 2 . d 



ti Jti dt J tl dt j tl 



dt — (vu)= j dt — (vu)+ / dt(— vu + vu) = vu . (2.5) 



ti 



ti 



There are two, linearly independent solutions to (2.3), which form the space of so- 
lutions associated with (2.3). In the case of constant positive A := u 2 , a basis for the 
2-dimensional solution space of (2.3) may be given by {sin(u;i), cos(ut)}. We shall remark 
here that, even for generic X(t), any two linearly independent solution curves behave like 
sine and cosine if they are to possess zeros. To this end, first we state the following lemma, 
which is easily verified by direct differentiation: 

Lemma 2. Let u(t) and v(t) be two linearly independent solutions of (2.3) and let it (to) = 
0. Then i>(to) 7^ 0, and u may be expressed as 



r ds 

u(t)=u(t )-v(t)v(t ) jf j^-^ (2.6) 

as long as v is non-zero. 

Formula (2.6) can be used to create a linearly independent solution with prescribed zero 
position and initial velocity to a given solution curve. Before embarking on further rela- 
tionships between v and u, we show 
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Lemma 3. Let v(t) be a (non-trivial) solution of (2.3) and let T\ be a time such that 
v(ti) = 0. Then either v(ti) ^ or else A has a pole at T\. In particular, if X is piecewise 
continuous, then the zeros of v are separated. 

Proof. Let to be a time where v(to) ^ 0, and define u(t) according to Lemma 2 with 
u(to) = 1. With the help of Lemma f applied to time to and t\ — e with e > 0, we have 

u(ri - e)u(ri - e) - v(t ) = u(ti - e)it(ri - e) . (2.7) 

Thus, by letting e go to zero, v(t\) = is possible only if u{j\) = oo or u{t\) = oo, 
which both imply that u can not be C 1 at T\. Therefore, A is neither continuous nor even 
piecewise continuous at ri, proving the assertion. Q.E.D. 

Let us henceforth assume that A be piecewise continuous or, if not, be at least such 
that whenever v{t\) =0, then v{t\) ^ 0. We then have the following 

Lemma 4. In the situation of Lemma 2 let T\ be the first zero position of v beginning 
with to. Then the formula (2.6) for u is well-defined at n and its value is given by 
u(ti) = -u{to)v{to) / v{t\) . 

Proof. Since u is a solution of (2.3) and A is piecewise continuous, any possible singularity 
of u stemming from its definition (2.6) can be restored, and the value of u at T\ can be 
calculated with the help of Lemma 1. One can check this also directly from (2.6) in the 
following way. Let eo and e be positive numbers, such that < e < eo < 1. Now if 
v{j\) = 0, then from (2.3) it follows v(j\) = 0. Thus in the interval [r\ — eo,T\] v can be 
approximated as v(t) = v(r\)(t — t\) + 0((t — ri) 3 ). Keeping eo fixed, but letting e go to 
zero we obtain 



u(n - e) = u(t ) (-eu(n) + C(e 3 )) v(t ) 



' r Ti - eo ds r i_£ ds 

J to R^ + X 1 - £0 (^ 



= 0(e) + u(f ) {-ev(n) + 0(e 3 )) ^o)j^^ £ dr |1 (l + 0{r 2 )) } 
= (9(e) - *fr°W* 0) , 
proving our assertion. Q.E.D. 

In the following proposition we now make precise in what sense two linearly indepen- 
dent solutions v and u behave like sine and cosine for the case when A is a positive constant 
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Figure 1. A schematic picture of how two linearly independent solutions, v and 
u, behave under a generic potential A. 

number. For the solution curves of the differential equation (2.3) to behave actually like 
sine or cosine, they must have zeros. This is certainly not the case for example for the 
constant A = — c 2 < with one possible solution curve being exp(ct). On the other hand, 
as we shall see further below, for the case of critical A, the solution curves necessarily have 
zeros. Thus for the following proposition we take it for granted that the solution curves 
indeed possess zeros, which are separated according to Lemma 3. 

Proposition 5. Let u and v be two linearly independent solutions of (2.3). Let t-\, 
to, ti, t 2 , ... be the zero positions of u, and r_i, tq, t\, t 2 , ■■■ be those of v. Assume 
to G (to, ti), v(t ) > and u(t ) > for definiteness (other cases of signs may be treated 
analogously). Then u is positive (negative) on (ti,ti + i) for i even (odd), v is positive 
(negative) on (r,, n+i) for i even (odd), and the zero positions of v and of u alternate such 
that 



Proof. First of all, the sign of u (and likewise of v) must alternate as one passes the zeros, 
since at each zero position the slope of u is non-zero by Lemma 3. We show first that 
the zero position ti following to nes m ( T i> T 2)- Since v is positive on (to,7~i), u may be 
defined by (2.6) and is therefore strictly positive on (t , T\). Moreover, since the derivative 
of v at the zero position n is non-zero (Lemma 3), and since v becomes negative when 
it crosses ri, we must have v(ti) < 0. Then, according to Lemma 4, u is strictly positive 
also at Ti, which implies T\ < t\. By applying Lemma 1 to the interval [ti,T2] one has 



... < T_i < t-i < Tq < t < Ti < ti < ... 



(2.8) 
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u{t2)v{t2) = u{t\)v{t\) = —u(to)v(to) < 0, the second equality following from Lemma 4. 
Since, by definition, v has a zero at r-i and becomes positive thereafter (until the next zero 
at T3), its derivative i)(j<i) is positive, which implies that ufa) < 0, i.e., u must have a 
zero at some t\ < t^- Having proven t\ <t\ < T2, we may proceed by re-defining u on the 
interval [ti,T2] analogously to (2.6), 

u(t) = u(t 1 )-v(t)v(t 1 ) jf . (2.9) 

The above reasoning for t\ may be applied again to yield T2 < £2 < T3 (see Fig.l). Pro- 
ceeding similarly for other ti, one may prove (2.8), Q.E.D. 

We remark that the formula (2.6) valid on [to, t{\ can be extended to the interval 
[Tj,Tj+i] by simply replacing the lower integration boundary to in (2-6) to U. To see this, 
one uses (2.9) and u{t\)v{ti) = u(t )v(to) (by Lemma 1). 

2.2. Critical potential and its characteristics 

We now consider the full equation of motion (2.1) on a fixed time interval [0, T], and 
examine the phenomenon of caustics. First we recall that, if u, v are the two linearly 
independent solutions of the homogeneous part (2.3) obeying the initial conditions, 

u(0) = 0, u(0) = l, u(0) = l, u(0) = 0, (2.10) 

then the general solution of the full equation (2.1) is given by 

x(t) = av(t)+/3u(t) + s(t) , (2.11) 

where s is a special solution of (2.1) which we we take to be one satisfying s(0) = 0. 
(We could further specify by requiring, say, s(0) = but this is not important at the 
moment.) The constants a, (3 are determined ^from the initial position x(0) and velocity 
x(0), respectively. Note that the solution u may be obtained from the Jacobi field (2.2) 
by u(t) = J(p,t)/j(p, 0). 

Caustics occur when the solution u(t), or the Jacobi field J(p, t), vanishes at t = T. If 
this happens then it follows that the final position x(T) = av(T) + s(T) does not depend 
on the initial velocity x(0) (which is also obvious from the definition of the Jacobi field). 
More precisely, we have 
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Lemma 6. Given a time interval [0, T], suppose that the Jacobi field J(p,t) beginning 
at t = vanishes at T. Then there exists a solution of the full equation of motion (2.1) 
satisfying the Dirichlet boundary condition, 

x(0) = a, x(T) = b (2.12) 

if and only if 

b = v(T) ■ a + s(T) , (2.13) 

where v(T) and s(T) are the final values of the aforementioned solutions. Further, for the 
final point b given in (2.13) there are infinitely many solutions satisfying (2.12), and each 
solution is uniquely characterized by the initial velocity x(0), which can take any value. 

Proof. In the general solution (2.11) the Dirichlet boundary condition (2.12) implies that 
a = a while (3 remains arbitrary since we now have u(T) = 0. Thus whatever the initial 
velocity (which could determine (3) may be, the final point x(T) turns out to be given by 
(2.13). Since (3 is left arbitrary, there are infinitely many solutions satisfying (2.12) with 
different initial velocities, x(0) = (3 + s(0). Q.E.D. 

Thus, if J{p, T) = 0, the different solution curves which spread out of the initial point 
a are all focused in one final point b. This phenomenon can be considered as a special case 
of caustics in geometric optics, see e.g., [5]. In this paper we fix the interval [0, T] and 
study the implications of caustics that arise under a generic potential A, and to this end 
we shall make the following 

Definition 7. If the Jacobi field J(p, t) vanishes at t = T, then the potential A shall be 
called critical, otherwise non-critical. 

In case of critical potentials, it is useful to introduce the constant 

MA) := ^ , (2 . 14 ) 

^from a solution v satisfying (2.3) with v(0) ^ 0, which gives the stretching factor during 
the period [0,T]. Clearly, this constant is independent of the choice of the solution v and 
hence determined solely by the critical potential A. Given an initial point x(0) = a, the final 
point (2.13) to which all solution curves are focused on now reads x(T) = k(X) a + s(T), 
and is called the conjugate point to a or focal point. The simplest example of caustics arises 
in the harmonic oscillator, where A = uj 2 is a positive constant at one of the frequencies 
to = ^ for n = 1, 2, — The stretching factor associated with the potential is then 
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k(X) = (— l) n , and it is easy to see that there are infinitely many solutions x(t) arriving at 
the conjugate point (— l) n a. 

The following Lemma shows how much the situation of a critical potential differs from 
that of non-critical one. 

Lemma 8. If A is non-critical, there exists a unique solution of (2.1) for the Dirichlet 
boundary condition (2.12) with any a and b. 

Proof. Since A is non-critical we have u(T) ^ 0. Then, for any a, b G IR, the Dirichlet 
boundary condition (2.12) determines uniquely the constants a, (3 in the general solution 
(2.11) to yield 

x(t) = a v(t) + b ~ aV ^ S{T) <t) + s(t) . (2.15) 

Q.E.D. 

We note in passing that a generic potential is non-critical, and this can be seen roughly 
as follows. Consider a non-critical potential and a solution curve of the force-free equation 
of motion (2.3) with initial value and initial velocity 1, which thus does not vanish at 
the final time T. If the potential is varied in an infinitesimal manner, then this solution 
curve also varies very little, implying that the final value remains nonzero for such small 
variations. This implies that the set of non-critical potentials is open in an appropriate 
topology. Further, by the same reasoning, if the potential is critical, then in general a 
small variation suffices to make the potential non-critical. Thus, again in an appropriate 
topology, which we do not specify here but is easy to find, the set of critical potentials is 
nowhere dense. 

The equation of motion (2.1) can be derived from the action functional I[x] = f Q dt L 
with the Lagrangian L in (1.1), which now acts on the space of paths x(t) with some given 
boundary condition x(0) = a, x(T) = b. If the potential is non-critical, then there exists 
a unique classical path from a to 6, denoted here by x(t), at which the action becomes 
a minimum. Suppose we vary the potential A and with it the classical path x(t) and its 
action I[x] in a small neighbourhood of a critical potential A with a and b fixed, wherein 
we only consider non-critical potentials. If A approaches A, two cases may occur: either b 
happens to be conjugate to a under A, i.e., b = ak(X) + s(T), in which case there exists a 
classical path also for the critical potential, and the action of course remains finite, or else 
b is not conjugate point to a, in which case there exists no classical path at the caustic, 
and the action becomes infinite, as the following lemma shows. 
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Lemma 9. In the situation of a non-critical potential A approaching a critical potential A, 
the minimum action I[x] of the classical path x(t) between two given fixed boundary points 
tends to infinity if the final point is not a conjugate one of the initial point for the critical 
potential A. 

Proof. Consider the classical path satisfying x(0) = a and x(T) = b given in the form (2.15), 
where we specify the initial velocity of the special solution s as s(0) = for definiteness. 
Denoting x(t) := av(t) + s(t), and using the differential equations satisfied by s and u, we 
obtain the identity, 

f T <M*) x(t) = f dt fiit) ~x{t) - (s(T)u(T) - s(T)u(T)) . (2.16) 

Jo Jo UK 1 ) 

With the help of this identity and Lemma 1 applied to u and v of (2.15), we find that the 
action becomes 

I[x] = ~\£ dtn{t) x(t) + l(2b- x(T)) x(T) + \ (b - x{T)f . (2.17) 

As A approaches A, we observe that u(T) — > 0, ii(T) — > l//c(A), x{T) — > k(X) a + s(T), 
while x(T) remains finite since its initial velocity x(0) = is fixed for all A. The assertion 
then follows ^from the assumption, b ^ k(X) a + s(T). Q.E.D. 

An important point to note is that the sign of the divergence I[x] — > ±oo depends on the 
sign of the product u(T) k(X) in the limit A — > A, and this depends on how A approaches 
A, not just on A it is approaching. 

If A is critical, then besides the stretching factor /c(A), a further characteristics is the 
index of the differential operator A\ viewed as a symmetric bilinear functional, which we 
now illustrate (see e.g., section 12 in [5]). Consider a generic potential A with external 
force fx and the action functional I[x] for paths x with fixed endpoints x(0) = a, x(T) = 6, 
and let x c \(t) be a classical path (which in case of critical A is assumed to exist) which 
obeys the boundary condition. To x c \(t) we add a perturbation n(t) with 1](0) = 1](T) = 
and expand the action I[x c \ + rj] in terms of ?y, the calculation of which is straightforward, 

I[x cl + n] = I[x cl ] + [ dtrj (A x x cl -fu) + l [ dt nA x r} . (2.18) 
Jo ^ Jo 

Since x c \ is a classical solution, the second term linear in r\ on the right hand side vanishes. 
The last term on the right describes the action of a bilinear functional on the space of 
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paths defined on [0, T] with vanishing boundary values. The index of this (non-degenerate) 
bilinear functional is then defined by the dimension of the space on which it is negative 
definite. It characterizes the type of the saddle point of the action functional I at the 
classical path x c \. 

From the consideration of Sturm-Liouville problem the operator A\ is known to pos- 
sess a complete set of orthonormal eigenfunctions u n (t) [13, 14]: 

u n {t) = E n u n {t) , (2.19) 
with 

u n (0)=u n (T)=0; [ dtu n {t)u m {t) = 8 nm . (2.20) 

Jo 

The index is then given by the number of negative eigenvalues of the eigenvalue problem 
(2.19) for non-degenerate A\ for which no zero mode exists. It is known that for large n 
the eigenfunctions u n approach trigonometric functions with the corresponding asymptotic 
eigenvalues (^-) 2 [13, 14], which are thus unbounded from above. In the following, we 
assume the eigenvalues to be bounded from below, which is for example the case when A 
is bounded from above. 

The bilinear functional A\ becomes degenerated if there arises a zero mode solution 
of (2.19), i.e., a u m with E m = for some m. The index of A\ may be extended even to 
this degenerate case by saying that it is given by the number of modes with E n < 0. By 
definition, the degeneracy occurs when the potential A is critical under the given interval 
[0,T]. To illustrate this fact from another viewpoint (see e.g., pp79 in [5]) let us consider 
the initial value problem x(0) = xq, x(0) = x$, where we no longer fix the final time T, 
but let it be variable. Looking for each T the classical path, the corresponding expansion 
(2.18), and the eigenvalue problem (2.19), with gradually growing T, we observe that the 
index of the operator A\ changes only at focal points when an eigenvalue goes through 
zero, i.e., at those values of T, where A restricted on the interval [0, T] becomes critical. 
The precise statement of this is the following 

Proposition 10. Let u be a solution of (2.3) with initial condition u(0) = 0. Then the 
index of the operator A \ is given by the number of zeros (focal points) of u on the half-open 
interval (0, T], which is called the Morse index. If X is critical and v a solution of (2.3) 
with the initial condition v(0) ^ 0, then v must have a zero on (0, T). The number of zeros 
of v is given by the Morse index. 



A\u n {t) = - 
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Proof. The first assertion on the index is a special case of the Morse theorem, cf. [15]. 
The second assertion follows directly from Proposition 5, since between two adjacent focal 
points of u there must lie a zero of v. Q.E.D. 



Note that, for A critical, the Morse index, which is characterized by A and is denoted 
by m(A), gives the number of negative eigenvalues of (2.19) plus one. We now have the 
following consequence of the foregoing discussions on the stretching factor of a critical 
potential. Let A be critical, and «bea zero mode solution of (2.3) with zeros at to = 0, t\, 
t m = T, where m = m(A) is the Morse index. Since a zero mode solution is unique up 
to a multiplicative factor, these zeros are fixed for a given potential A. Let Aj := A|[ ti _ ljti ] 
be the restriction of A on the i-th interval [£;_i,£i], i = 1, m. Then the restriction of 
the zero mode u on each such interval is trivially a zero mode solution corresponding to 
Aj. Thus Aj is critical with Morse index m(Aj) = 1. Then one has 

Lemma 11. The stretching factor of a critical potential is the product of the stretching 
factors of the potentials restricted to each closed interval between two adjacent zeros of a 
zero mode solution, i.e., 



Each k(Xi) is negative, and thus the stretching factor is negative (positive) if the Morse 
index m(A) is odd (even). 

Proof. Let v be a solution of (2.3) linearly independent from it. Then from the proof of 
Lemma 6 we have k(Xi) = v{U) /v(U-i) and 



Each stretching factor is negative, since by the proof of Proposition 5 between two adjacent 
zeros of the zero mode solution u there must lie a zero of v, where its sign changes. Q.E.D. 



k(X) = k(X 1 )---k(X m ). 



(2.21) 



k(X) 




(2.22) 
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3. Quantum Caustics 



In the previous section we studied various dynamical aspects of a particle moving 
under the influence of time dependent harmonic potentials X(t) and external driving forces 
/j,(t). When X(t) is critical, the dynamics of the system exhibits a number of singular 
characteristics at the classical level, notably in that transitions are allowed only between 
conjugate points specified by the stretching factor k(X) of the potential. The aim of this 
section is to study the quantum dynamics of the system by looking at the transition 
amplitude closely. In particular, for critical A(t), we shall confirm the known fact [5, 6] 
that classically forbidden transitions continue to be forbidden even at the quantum level. 
The novelty of our result is the closed form of the path-integral kernel for the transition 
amplitude expressed in terms of the stretching factor k(X) and the Morse index m(A). We 
further discuss the quantum effect in the momentum susceptibility of a Gaussian wave 
packet for non-critical X(t). 

3.1. Transition amplitude and caustics 

We begin by considering the transition amplitude for the system defined by the 
quadratic Lagrangian (1.1). Let H(t) be the Hamiltonian operator corresponding to the La- 
grangian, and U (T, 0) the unitary operator for the time evolution obeying the Schrodinger 
equation, (id/dt)U(t,0) = H(t)U(t,0). Then, in the path-integral formalism, the ampli- 
tude for the transition from x = aa,tt = Otox = ba,tt = T, where a and b are arbitrary 
two points on the line (R, is given by 



In carrying out the path-integration, we need to take into account the fact that, when 
X(t) is critical, there may not exist a classical solution that respects the given boundary 
condition. The general case, including the critical one, may be handled by the following 
procedure [5, 4]. 

First, let c be the endpoint at t = T of a classical solution x c \ of the equation of 
motion (2.1) starting from a at t = 0, i.e., x c \(0) = a and x c \(T) = c. Here we put the 
bar to the solution in order to emphasize the fact that the endpoint value c may not be 
equal to 6, reserving the notation x c \ without a bar for the actual solution, if any, having 
x c \(T) = b. If A is non-critical, Lemma 8 states that c can be chosen arbitrarily and hence 
the solution x c \(t) does exist, whereas if X(t) is critical, Lemma 6 states that the endpoint 




(3.1) 
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is determined uniquely c = k(X) a + s(T) by the stretching factor of the potential and the 
endpoint of the solution s(t). Now, let us choose a fixed, smooth function p(t) satisfying 



p(0) = 0, p(T) =b-c. 



(3.2) 



Then we may decompose any path x(t) connecting the endpoints, x(0) = a and x(T) = 6, 
as 

x(t)=x cl (t) + p(t)+ V (t). (3.3) 

The function 77(f) in (3.3), which represents the fluctuations, vanishes at both of the ends 
7/(0) = rj(T) = and may be expanded rj(t) = J2 n a n u n (t) in terms of the orthonormal 
eigenfunctions in (2.19) and (2.20). 

Using the decomposition (3.3) we find 

I[xd + P + V]= I[xd] + I[p] + I[v] 

+ / dt{x c ip + pf] + r]x c i - X(x c ip + pr] + r]x c i)} 
Jo 

= I[x cl ]+I[p}\»= + lj2 E " a n (3 ' 4) 

n 

+ x c i(T)p(T) + p(T)^2a n u n (T) + ^2E n a n / dt p(t)u n (t) . 

n n J° 

After the change of the integral variables T>x = T>rj = M \\ da n where M is a Jacobian 
factor, the path-integral (3.1) becomes a product of Gaussian integrals over the modes a n . 
If A is non-critical, there exists the solution x c \ and hence we are allowed to choose p(t) = 
identically in (3.4). The Gaussian integrations then gives the standard result, 



K(b, T; a, 0) = A/" 



(3.5) 



which can be shown [2, 3] to be equivalent to (1.2). 

If, on the other hand, X(t) is critical, then by definition we have E m = for some m. It 
follows that the integration over a m , which is now non-Gaussian, yields the delta-function, 
S(p(T) u m (T)). Hence, for critical potentials the kernel formula (3.5) should be modified 
into 

'2tt 



K(b,T;a,0) 




\{E n 



5(p(T)u m (T))e u ^\ 



(3.6) 
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where we have set p(t) = in the phase part of the kernel (3.6), which is allowed in the 
presence of the delta- function and due to the fact u m (T) ^ 0. Since we have now 



p(T) = b - k(X)a - s(T) , 



(3.7) 



the kernel (3.6) already shows that classically forbidden transitions to non-conjugate points, 
b 7^ k(X)a + s(T) from a, are forbidden even quantum mechanically. We remark at this 
point that, strictly speaking, the action I[x c \] in (3.6) should be written as I[x c \] since 
there is no solution x c \(t) if p(T) ^ 0. However, the present form is still possible because 
the phase becomes insignificant if p(T) ^ in which case the kernel vanishes identically. 

The remarkable fact in the critical case is that one can express the transition kernel 
(3.6) in terms of the stretching factor k(X) and the Morse index m(A) of the potential X(t). 
To see this, let us first write the kernel (3.6) in the form, 

K(b,T;a,0) = R(T)S(p(T))e i ^ T \ (3.8) 

with 

R(T) = y/2^\Af 

Then, the unitarity relation, 

J dbK*(b,T;c,0)K(b,T;a,0) = J db (c\U^(T,0) \b)(b\U(T,0) \a) = 6 (a - c), (3.10) 



n W 

n^=m 



\Um(T))\ 



-1 



(3.9) 



determines the modulus R(T) to be 



R(T) = ^\k(Xj\. 



(3.11) 



On the other hand, since the Morse index gives the number of negative modes plus one, 
we have 





2 




\[E n 




n i^i 









-f«A)-l) 



e 2 



Thus the phase part is given by 



TV 



Q(T) = I[x cl ] - -(m(X) - 1) . 



(3.12) 



(3.13) 



By combining these we find that the transition kernel (3.6) takes the simple form 4 

K(b, T; a, 0) = v^(A^5(6 - k(X)a - s(T)) e iiM-^m(X)^ (3J4) 



An abstract form of the formula for the period before the first conjugate point has appeared in [6]. 
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up to an overall constant of unit modulus. 

As application, take, for example, the forced harmonic oscillator given by X(t) = u 2 
and ji{t) = —f(t). A special solution satisfying s(0) = is then 

1 f* 

s (t) = - dt' smoj(t -t')f(t '). (3.15) 
w Jo 

As mentioned earlier for the pure harmonic oscillator case, caustics occur at uj = ^ with 
n = 1,2,..., where we have the stretching factor k(X) = (— l) n and the Morse index 

rp 

m(A) = n. The endpoint of the solution is s(T) = — (— l) n f Q dt smut f(t)/oj and the 
classical action reads 

I[x cl ] = a [ dt cos utf(t) - - f dt [ dt' cosut smut' f(t) f(t') . (3.16) 
Jo v Jo Jo 

In particular, for constant / the kernel (3.14) reduces to 

K(b, T; a, 0) = s(b - (-l) n a - {1 - (-l) n }//w 2 ) e ^l^~^m, (3.17) 

This agrees with the results obtained earlier [16, 17, 18] by other indirect means 5 for the 
(forced) harmonic oscillator. 

It is possible to combine the two kernel formulae, (3.5) and (3.6), into one expression. 
For this, we just introduce an infinitesimal e > to replace E n with E n + it in (3.4), and 
perform the integrations over a n (which are now all Gaussian for any A) to get 



K(b,T; a,0) = M 



Y[(E n + ie) 



e i*(6,a;A) ^ (3 ^ 



where 

$(6, a; A) := 7[x c i] + I[p]\^o + ^i(T)p(T) 



(3.19) 



Now for A non-critical, the kernel (3.18) reduces to (3.5) in the limit e — > by construction. 
(This can also be confirmed explicitly; see Appendix.) But it is also easy to see that, for A 
critical, it reduces to (3.6) in the limit, since the mode m with E m = yields the required 
delta-function on account of the identity, lim e ^o 77^ e - x2 / 2e = S(x) . This expression 



The result contained in [18] for f(t) 7^ is marred by an error in the phase. 
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(3.18) with (3.19) will be useful later in section 4 when we examine the correspondence 
hypothesis. 



3.2. Quantum effect in the Gaussian slit experiment 



For A non-critical but close to a critical A, one expects that there will be a concen- 



tration in the transition amplitude around the focal point b = k(X)a — s(T) conjugate 
to a given initial point a. We now analyze how the concentration takes place quantum 
mechanically, based on the Gaussian slit experiment which is laid out in the textbook of 
Feynman-Hibbs [19]. 

To illustrate our point, let us first consider the harmonic oscillator A = uj 2 (with 
[i = 0). In the Gaussian slit experiment we look at the evolution of a Gaussian wave 
packet prepared at t = 0, and for this we put an apparatus which emits a particle from the 
origin x = at time t = — r. To get a Gaussian distribution at t = 0, we place a 'Gaussian 
slit' centered at x = a with effective width (variance) oo- The slit is realized by means of 
the Gaussian transmission factor, 



The wave function ifj(x,0) of the particle at t = is then furnished by the product of the 
free particle kernel 6 K (x, 0; 0, — r) = ^/l/(2nitiT) e lx2 / 2hT and the transmission factor. 




(3.20) 



From the normalization condition J dx \ip(x, 0) | 2 = 1 we determine the constant iV and 
obtain 




Note that the initial state (3.21) has the average momentum 




(3.22) 



At t = T, the wave function is given by 




(3.23) 



6 



To render the quantum effect manifest, we keep ti throughout this subsection. 
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with the kernel for the harmonic oscillator [16], 



2nih\ sinu/Tl 



exp 



2ft sin 



{(x 2 + y 2 )cosuT-2xy}) e"^ Int [ 



K(y,T;x,0) 

(3.24) 

where Int[x] stands for the greatest integer below x. Together with (3.21) the integration 
in (3.23) is Gaussian and hence we get 



(27Ri 2 (T))i 



• , (// , ' l ':; i i n] ■ / (phase) 



4<x 2 (T) 



(3.25) 



Thus, as is well-known, the Gaussian distribution retains its shape for the harmonic os- 
cillator, in which the center moves along the classical path having the initial momentum 

P, 



x c \ (T) = a cos ujT H sin uT, 

UJ 



(3.26) 



while the variance pulsates as 



i. 



a 



(3.27) 



Clearly, the second term in (3.27) represents the quantum effect whereas the first term is 
just the classical variance, since the variation in the initial position Aa; c i(0) = o~o around the 
center x c i(0) = a for the classical path results in the final variation Ax c i(T) = (7q\x c \{T)\/ a. 
(Notice that the initial momentum p is linearly dependent on the initial position where the 
path goes through.) Thus, as one expects, the quantum effect always enhances the spread 
of the Gaussian distribution. 

Viewed as a function of the initial variance ctq, the final variance cr(T) in (3.27) attains 
its minimum, i.e., the highest concentration of intensity, 



at 



hx cl (T) sinwT 



ha sinwT 



(3.28) 



(3.29) 



2uox cl (T) 

which is precisely the point where the quantum effect matches the classical variance. The 
infinite concentration a m \ n {T) = takes place at u> = nrr/T, or at p = —au/ tancuT. The 
former case is the expected caustics, where the variance (3.27) becomes purely classical 
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and reproduces the initial value c(T) = cr . Thus the infinite concentration is obtained 
if we let (To — > 0. By contrast, in the latter case the variance becomes purely quantum 
mechanical, as all paths passing through the slit coalesce toward the origin x = at t = T. 
In fact, this case is again caustics occurring during the combined period [— r, T] under the 
potential \{t) which vanishes for t < 0. Note that in this case the infinite concentration 
at the quantum level is achieved by letting o~o — > oo. In practice, however, both of these 
concentrations in intensity are unstable, since a small fluctuation in the parameters u> or 
t will generally bring the variance to a large value, as can be seen in Fig. 2. Nevertheless, 
one could achieve a very high intensity by adjusting the initial variance along with the 
parameters according to (3.29). 

An important quantity for characterizing the concentration is the susceptibility of the 
variance against initial momentum fluctuations. Exposing the momentum dependence of 
the variance explicit <r(T) = a(jp,T), we use the following normalized quantity for the 
susceptibility, 

S(p f T) := ~a(p t T). (3.30) 
For the present harmonic oscillator (3.27) we find 



S( P , T) 



sin uT 


H 




UJ 





UoWtJ • (3 - 31) 



In the classical limit % — > 0, the susceptibility reduces to the (absolute value of the) 
Jacobi field, which in this case becomes J(p,T) = smujT/uj. In the two cases mentioned 
above in which an infinite concentration can in principle be possible, the susceptibility 
vanishes S(p, T) = and hence the Gaussian wave packet becomes free from momentum 
fluctuations. What is interesting in this result (3.31), however, is that the quantum effect 
suppresses the susceptibility of the variation against initial momentum fluctuations. 

We now show that these features are universal for a generic potential A if \x = 0. For 
this we need to note two crucial points. First, for systems with quadratic Lagrangians the 
semiclassical (WKB) approximation for the kernel gives the exact result, 



K{y,T;x,0) = 



2nh 



d 2 L 



cl 



dydx 



2 



,'J cI -^m(A) ^ (332) 



where I c \ = I c i(y,T; x,0) := I[x c \] is the action evaluated for the classical path x c \(t) 
connecting the endpoints, x c i(0) = x and x c \(T) = y. 
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Fig. 2b 




G( c ,p,T) 



Figure 2. The final variance <r(T) = a(ao,p, T) as a function of the initial variance 
(To and the momentum p. Fig. 2a shows a generic case given by the parameters 
a = 0.3, u = 2.5 and T = 1.0 (we set h = 0.004 here). The focal points due to the 
caustics associated with the period [— r, T] appear at the bottom of the valley at 
p = 1.0. The minimum value a m i n (T) at the bottom tends to vanish as o~o —> oo 
but it is unstable under the fluctuation of p. When u) reaches the first critical value 
u = it associated with the period [0, T], the situation changes drastically to the 
one shown in Fig. 2b. This time the minimum o" m ; n (T) = attained at do = is 
unstable under the fluctuation of u. 

Second, if \x = the classical action I c \(y,T;x,0) is a quadratic polynomial homoge- 
neous in x and y, because we have 

I c \{cy, T; cx, 0) = I[cx c i) = c 2 /[x cl ] = c 2 / c i(y, T; x, 0) , (3.33) 

for any constant c. Thus, with A, B and C being some functions of T, we may write 

I c i(y,T;x,0) = Ax 2 + Bxy + Cy 2 . (3.34) 
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Then it is straightforward to see that the Gaussian wave packet (3.21) evolves into the 
wave function (3.23) at time t = T in exactly the same form (3.25) with 



x cl (T) = --(2aA + p), 



(3.35) 



and 



The variance attains its minimum, 

hx c i(T) 




(3.36) 



Cmin(r) 



aB 



at ctq 



ha 



2Bx cl (T) 



(3.37) 



The susceptibility of the variance (3.30) then reads 



S(p,T) = \J(p,T)\{l + 



ha 



2a$Bx cl {T) 



(3.38) 



where J(p,T) = — (d 2 I c \/dydx) 1 = —1/B is the Jacobi field for the classical action 
(3.34). The generic behavior of the susceptibility (3.38) with respect to p is analogous to 
what we have seen in the harmonic oscillator case. 

The above features of the Gaussian wave packet persist even for fx ^ 0, because 
then the classical action I c \(y, T; x, 0) acquires only linear and constant terms in addition 
to the quadratic terms in (3.34) as can be explicitly seen from the form of the general 
classical solution (2.15). Accordingly, the time evolution by the integral (3.23) is essentially 
unchanged. Semiclassically, this may also be the case for more general systems, not only 
for those with quadratic Lagrangians we considered, in view of the earlier study [20] which 
suggests that these features are a norm for a generic system in the limit h — > 0. 
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4. The Role of Instantons in the Quantum Mechanical Rotor 



In the previous section we obtained the path-integral kernel for the quadratic action 
with a generic time-dependent potential function X(t). The aim of this section is to apply 
it to the problem of quantum mechanical rotor on the circle S 1 and thereby examine 
Jevicki's correspondence hypothesis [12] mentioned in the Introduction. We first furnish 
a general scheme to compute the path-integral kernel for the rotor based on the formula 
we just obtained. This will strengthen the basis on which Jevicki's arguments for the 
hypothesis stand and highlight their loose ends at the same time. Our close examination 
will show that neither of the two crucial assumptions adopted in his arguments can be 
entirely justified for generic transition amplitudes. 

4.1. Examination of the correspondence hypothesis 

The system taken up by Jevicki to illustrate the correspondence (1.3) between sta- 
tionary points and simple poles is the quantum mechanical rotor, namely, a free particle on 
the circle S . Embedded in the two dimensional plane IR , the system acquires a harmonic 
potential and develops a singularity in the transition kernel, and this has been used as 
a prototype to argue that the presence of instantons can be seen as simple poles (rather 
than stationary points) in the <LP N nonlinear sigma model. A salient feature of the S 1 
system is that its non-trivial topology — S 1 is multiply-connected: 7r 1 (5' 1 ) = Z — admits 
multi-winding configurations called 'instanton' solutions labelled by the winding number 
n G 7L. Here, our concern lies in the question whether the instantons can also be seen 

1 2 

explicitly even if we embed the S in the topologically trivial plane IR using the Cartesian 
coordinates. Before enforcing this 'Cartesian point of view', we wish to recapitulate the 
known result of the path-integral on S 1 in the 'intrinsic point of view', that is, when the 
circle is regarded as intrinsic rather than being embedded in (R 2 , where the instantons arise 
as sole contributors to the transition kernel. 

For a free particle on 5" 1 coordinated by the angle <p G [0, 2tt) we have the classical 
action 1$ [tp] = dt^(ip) 2 . In the path-integral formalism, the transition amplitude going 
from a at t = to f3 at t = T where a, (3 G [0, 2tc) is given by 



*/<p(0)=a 

This formal expression must be given a meaning such that it accommodates multi- winding 
paths allowed for the transition since S 1 is multiply-connected. The conventional method 




(4.1) 
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[5] to deal with such paths is that, instead of working on S 1 , one considers the corresponding 
free particle system on the covering space of the circle, i.e., the line IR governed by the 
action Iq[x] = J Q T dt^(x) 2 , which is invariant under the translation by 2tt. Namely, to the 
amplitude on S 1 consisting of paths of winding number n, one assigns the free particle 
amplitude K (f3 + 2nn : T; a, 0) on the line IR, and then sums it up over the integers n£2 
with a weight factor w(n), 

oo 

Kf(0,T;a,O)= w(n) K (/3 + 2nn,T; a,0). (4.2) 

n= — oo 

In short, this amounts to regarding 5" 1 as the coset space IR/Z by identifying those 
points on the line (R which differ by 2tt x integer. The weight factor w(n) appearing in (4.2) 
signals the ambiguity that can arise due to the multiply-connectedness of the space. It is 
furnished by the unitary representation w(n) = e in9 of the group Z, where 9 G [0, 2ir) is the 
angle parameter specifying the representation and, hence, the ambiguity. In what follows, 
however, we put 9 = for simplicity. Recall that on the line IR the kernel for the free particle 
is K (b,T; a, 0) = j^^Vxe 1 ^ = y/l/(2iriT) e u °^ with/ [x c i] = (b-a) 2 /2T. In our 
case, the classical solution satisfying the required boundary conditions on 5 1 possessing 
winding number n, the n-instanton solution, is given by 

^ (t) = ((3-a + 2n7i) | + a . (4.3) 

Accordingly, the kernel (4.2) on the circle 5" 1 becomes 



Kt 03, T; a, 0) = J -±- £ e ^~^ f /2T (zU) 



2niT 

n=—oo 



which is just a sum of instanton contributions. 

In the Cartesian coordinate point of view, on the other hand, one puts the particle 
on the two dimensional plane coordinated by x = (x, y) G IR 2 and imposes the constraint 
|x| 2 = 1. The constraint can be implemented by introducing the Lagrange multiplier \{t) 
and writing (4.1) as 

i / .x(T)=(cos/3,sin/3) • rT 

Kg (p,T;a,0)= VXVxexp- dt (|x| 2 - A(|x| 2 - 1)) 

Jx(0) = (cos a,sina) ^ JO ^ 



I 



i r T MT)=cobP f y(T)=sin(3 

VXe^Jo dtx Vxe lI[x] I Vye lI[y] 

Jx(0)=cosa ^y(0)=sina 



(4.5) 



For a general discussion of the path-integral quantization on a coset space G/H, see [21]. 
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where I[x] (and similarly I[y\) is the quadratic action with time-dependent harmonic po- 
tential \{t) under zero force fi(t) = 0. In words, the kernel is given by the product of two 
kernels on the line sharing the potential A(£), one of which describing the transition ^from 
cos a to cos/3 and the other from sin a to sin/3, averaged over X(t) with the weight factor 

/■T 
dtX 
o , 

Kf(p,T;a,0) = J VX K(cos/3, T; cosa, 0) K(sin0, T; sina, 0) e* -A dtX . (4.6) 
Using the expression (3.18) for K{b, T; a, 0) obtained in section 3, the kernel (4.6) becomes 

r -1-1 

Kf(p,T;a,0)=fvX ]\(E n +ie) e "Kcos/3,cosa;A)-M^sin,Mna;A) + i dt A _ ^ 

L n 

Now the generic formula for the phase (3.19) implies that the phase $(&, a; A) reduces 
simply to the action I c \(b, T;a,0) = I[x c \] except for the cases where the potential A is 
critical and b is not conjugate to a, i.e., b ^ k(X)a. However, these exceptional cases may 
be neglected in the integration over A on the grounds that they form a set of measure zero 
in the entire space of potentials A and that the kernel K{b, T; a, 0) vanishes at those A and 
hence cannot contribute to the integral (4.6). For this reason we replace the phase factor 
in the integration by the classical action with the prescription, 

$(6, a; A) - 7 c ( f } (6, T; a, 0) := (1 + i sign(J)£) I cl (b, T; a, 0) . (4.8) 

The prefactor (1 + i sign(/)£), where sign(J) := I/\I\ gives the sign of the action with £ > 
being infinitesimal, is attached in order to guarantee that e^ 1+isign ^^^ vanishes for those 
critical potentials with b ^ k(X)a for which the action diverges (see Lemma 9). 

To proceed further and make contact with Jevicki's argument, let us consider the 
separation of the zero mode part Ao ^from the potential, 

X(t) = X'(t) + A , with I dtX'{t) = 0, A = ^ / dtX(t) . (4.9) 

Jo T J 

As we did for the potential A(£) in (2.19) and (2.20), we can also consider a complete set of 
orthonormal eigenfunctions for the non-zero mode part X'(t). Let the eigenvalues of those 
eigenfunctions be e n . Obviously, both X(t) and X'(t) share the same set of eigenfunctions 
with eigenvalues related by 

E n =e n - A . (4.10) 
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Accordingly, a potential X(t) is critical if A = s n (X') for some n. 
Substituting (4.10) in (4.7), we find 

r n _1 

Kf (/?, T; a, 0) = / V\'d\ \[{e n - A + ie) 

^ In 

x exp ji^P (cos P, T; cos a, 0)+z/ c ( P(sin/?, T; sin a, 0) + -^qT 

(4.11) 

Consider now the integration over the zero mode Ao- Observe that the integrand diverges 
(in the limit e — > 0) for critical A, because of the factor [n n ( £ ™ — + ^ e )] _1 - Recall 
that the additional divergence that could arise in the actions 7 c i(cos/3, T; cos a, 0) and 
2ci(sm T; sinct, 0) is taken care of by the prefactor + i sign(/)£). Thus we are naturally 
led to evaluate the integration over the zero mode by collecting the values at the critical 
potentials which admit the endpoint b to be conjugate to a in both x and y-direction. This 
will be carried out by regarding the Ao integration as part of the contour integration in 
the complex A plane as shown in Fig. 3. Since each critical potential which satisfies the 
above condition provides a simple pole in the integrand in (4.7), and since the contribution 
^from the semicircle of the contour vanishes as the radius becomes large due to the weight 
factor e^ A ° T , the Ao integration boils down to the residue calculus at those simple poles. 
That the Ao integration in (4.7) can be carried out by the contour integration in the way 
described above, ignoring the contributions from those critical potentials which do not 
meet the condition b = k(X)a, constitutes our first assumption. 

Unfortunately, this assumption cannot be true in general. Indeed, the condition for a 
critical potential X(t) to create just a simple pole (not other type of potential singularities) 
is to have the stretching factor k(X) which matches the boundary conditions in both x and 
y-direction simultaneously, 

, . , . cos 8 , , _ . sin 8 , t „ _ 

HA) = and HA) = , 4.12 

cos a sin a 

but, obviously, these are incompatible in general. It then follows /from the absence of 
simple poles that the transition amplitude (4.7) on the circle must vanish for generic 
boundary conditions, a result in contradistinction with (4.4). Having no simple poles 
that could correspond to the stationary points in the intrinsic viewpoint, we find that the 
correspondence hypothesis (1.3) cannot be true in general. It is also worth mentioning 
that in the intrinsic viewpoint the stationary points contributing to (4.4) do not in general 
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Im X 




Figure 3. The integration contour in the complex Ao-plane. Among the potentially 
singular points Xq = e n + it, those critical potentials which meet the condition 
stated in the text give rise to simple poles shown by By contrast, those which 
do not meet the condition, shown here by 'x', will not create actual singularities 
for contributing to the integration. 

correspond to the stationary points appearing under the critical potentials in the Cartesian 
viewpoint. 

The condition (4.12) becomes compatible if tana = tan /?, that is, if the two endpoints 
are the same or diametrically opposite with each other. To examine whether it is possible 
to sustain the correspondence hypothesis if restricted to this specific case, let us take 
a = f3 = from now on. Then we observe that, in this case, we have / c i(0, T; 0, 0) = for 
the classical action in y-direction, and the condition is fulfilled if k(X) = 1. We thus arrive 
at 

Kf (0, T; 0, 0) = F(T) f VX' e**" (V) , (4.13) 

J n 

with the total phase 

(4.14) 



*n(V) := (l^\l,T;l,0)-m(X)n + ^X T 



Aq— e r 



where we formally included all the contributions from critical potentials using (3.14) in 
view of the fact that those which correspond to A with k(X) ^ 1 will drop out due to 
the prefactor attached in (4.8). In (4.13), the real function F(T) accounts for the possible 
overall factor that arises in the change of integration measure and the residue computation, 
and m(A) is the index of the critical potential X(t). An important point to note is that, 
for any critical potential with k(X) = 1 that can contribute to the kernel (4.13), the index 
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to (A) is always even on account of the alternate nature of the sign of the stretching factors 
k(X) (see Lemma 11). Thus we may drop the index term m(A)7r in the total phase $ n (A') 
without affecting the kernel (4.13). 

According to Jevicki [12], we now pose the second assumption: the total phase factor 
<J? n (A') in (4.13) is independent of X'(t). This is a rather drastic assumption, allowing us 
to set X'(t) = and getting the eigenfunctions and the eigenvalues, 

/ \ /~2~ . / mvt\ /mv\ 2 . A „ „. 

un(t) = y-sin^— J , £ n={-Y) , (4-15) 

for n = 1, 2, 3, . . .. We then find that the critical potentials arise at X(t) = Xq = e n , and 
they have the stretching factor k(X = e n ) = (— l) n . Thus, among them those possessing 
k(X) = 1 occur at even integers n = 21. Since we have J c i(l,T; 1, 0) = for the classical 
solution in the (critical) harmonic oscillator, and since the index is given by to(A = 621) = 
2/, the total phase becomes 

$ 2Z (0) = -2/7T+^^ . (4.16) 

Consequently, we obtain 

Kl 1 (0, T; 0, 0) = -F(T) ^ ' T . (4.17) 

Thus, if we choose F(T) = 2^/^^, the result (4.17) coincides with (4.4), up to the 
contribution from / = 0. 

The foregoing argument is (a rigorous version of) the one used in [12] which asserts 
that, even in the Cartesian point of view, one can see the instantons contributing to 
the kernel by the form of simple poles, rather than stationary points, in the integrand. 
In fact, for the boundary condition a = (3 = the classical motions allowed under the 
critical potentials sitting at the simple poles correspond precisely to the instanton solutions 
(4.3) which are the stationary points in the intrinsic viewpoint. Thus we find that the 
correspondence hypothesis (1.3) holds if one restricts oneself to the specific transition 
process and ignores the lacking piece / = in (4.17) which, if existed, would correspond to 
the trivial (0-instanton) solution x^\t) = in the intrinsic viewpoint. 

4.2. Validity of the second assumption: an example 

Although we have seen that the first assumption cannot be sustained for generic 
transition processes, at the end of the computation the kernel (4.7) must recover the 
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integers (i.e., instanton numbers) appearing in (4.4) for any a and (3. It is a matter 
of fact that the classical motions expressed by the instantons x^\t) in (4.3) arise at 
X(t) = constant, even though they do not necessarily appear as simple poles for generic a 
and p. In this respect the second assumption does lead to the correct result in effect, and 
one is curious whether there is an a priori reason for this. (Note that the two assumptions 
are independent of each other.) In the rest of this section we shall investigate this issue by 
looking at a simple but nontrivial example. 

To this end, we return to the special case a = f3 = where the assertion of the 
first assumption is valid (up to the contribution at / = 0). The second assumption now 
turns out to be that, modulo 2-7T, the total phase $ n (A') in (4.14) is given by the phase 
&n(X' = 0) in (4.16) with n = 21. To make the issue more transparent, let us fix some \'{t) 
and consider the family of potentials \{t) which differ by constants X(t) = X'(t) + Ao as in 
the form (4.9). Within this family, the critical potentials which arise at A = e n for n G 7L 
form a class and can be labelled by the integers n. In particular, those that have k(X) = 1 
form a subclass in the class and may further be labelled by another set of integers / G Z. 
The second assumption can then be justified if the total phase $ n admits the form (4.16) 
modulo 2% with the summation / of n = 21 being performed over the critical potentials in 
the subclass. 

To see how this works, take the potential 

X'(t)=c5(t-t )-^ , (4.18) 

where c is a constant and the term — ^ is inserted to ensure J dtX' = 0. This is a pertur- 
bation of the trivial case X'(t) = in the sense that (4.18) defines a one-parameter family 
of potentials containing the trivial case c = 0. The other parameter to G (0, T) specifies 
the moment of impact inflicted on the particle and may be determined by requiring that 
X(t) admits cases where k(X) = 1. From Lemma 6 we know that this requirement is met if 
there exists a classical solution u(t) of the equation of motion vanishing at the endpoints 
it(0) = u(T) = and further has the same velocity -u(0) = u(T). Below we shall seek 
for conditions under which such a solution exists, setting the initial velocity ii(0) = 1 for 
definiteness. 

For this purpose, we first divide the interval [0, T] into [0, to] an d [to,T], and call 
them interval (I) and (II). Then, in each of the intervals, the solution which respects the 
boundary conditions takes, respectively, the form 

1 1 

m(t) = - 8m(ut) , u u (t) = sin(w(T - t)) . (4.19) 
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The continuity condition tti(to) = ^ii(to) at t = to gives 

1 1 

- sin(wto) = sin(w(T - t )) . (4.20) 

The velocity acquires a gap at t = to prescribed by the delta-function interaction as 
^n(to) _ ^i(to) + c ui(to) = 0, which amounts to 

cos(cu(T - t )) - cos(cuto) + — sin(u;to) = . (4.21) 

The first condition (4.20) is fulfilled if 

U (T - t ) = -vt + 2irl , or u(T - t ) = ut + tt(2Z + 1) , (4.22) 

for I G Z. To each of these cases, the other condition (4.21) requires that (for c ^ 0) 

2tvI m 
u = — , c = arbitrary, t = — T, (4.23) 

where m can be any integer with 1 < m < 21 — 1, or 

(21 + 1)tt 2ujcos(ujto) „„,, 

^ = L » c = wTT » to = arbitrary . 4.24 

T — 2t sm(o;to) 

Note that the conditions (4.24) do not allow for a set of infinite number of solutions for 
fixed c and to (specified by A'(t) that defines the class). Hence in the latter case (4.24) the 
subclass with k(X) = 1 consists of finite number (possibly one) of critical potentials. In 
the former case (4.23), on the other hand, if m = with p and q being coprime positive 
integers satisfying ^ < 2, then the parameter t lies in the interval < t < T with the 
fixed value, 

to = ±T. (4.25) 

Thus the former case does provide the subclass consisting of infinite number of potentials, 
but for m to be integers the label / runs as I = p, 2p, 3p, . . .. We now examine whether 
our second assumption holds or not here. 

To evaluate the total phase explicitly, we consider the classical solution x(t) for the 
above class of potentials satisfying the boundary condition, 

x(0) = x(T) = 1 , x(0) = 1 . (4.26) 
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(Again, for definiteness we added the second condition in (4.26) because the action on 
caustics does not depend on the initial velocity.) As before, in each of the intervals (I) and 
(II) the solution obeying the boundary conditions (4.26) takes the form, 

xi(t) = — sin (art) + cos(a;t) 

u (4.27) 
XlI (t) = -- sin(w(T - t)) + cos(w(T - £)) , 

where 7 is a constant. ^From the continuity condition and the velocity gap equation at 
t = t in (4.25) together with the 00 in (4.23), it is determined to be 7 = 1 — c. For this 
solution (4.27) the classical action becomes 

7 cl (l,T;l,0) = ^xi(t)xi(t)|* + ^xu(t)x n (t)\ T to - ^xj(t ) = ~ • (4.28) 
On the other hand, the eigenvalues for the critical potentials in the subclass are 

2tt/\ 2 c 



£2l= \~Y) + T' (429) 

where we put the label 21 to match the notation (4.15). Consequently, the total phase 
(4.14) turns out to be 

$ 2Z = -m(X)7r + . (4.30) 

Note that the final expression (4.30) is independent of the parameter c. Since the index 
m(A) is always even as remarked earlier, we see that, if p = 1 for which the label / takes 
all integers, the phase factor agrees with the original one (4.16) even if it is perturbed by 
the delta-function interaction. However, for p > 1, the total phase does not reproduce the 
unperturbed value (4.16). 

We thus find that the assertion of the second assumption holds for the family of 
potentials furnished by (4.18) with the parameters given by (4.23), if one makes the special 
choice p = 1. It does not hold, however, for the case (4.24), where the final phase factor 
consists of a finite number of contributions. In conclusion, there is no a priori reason that 
the second assumption can lead to the correct transition amplitude, simply because there 
are cases where the assumption itself breaks down. Finally, we point out that the case 
p = 1 where the assumption does hold implies to = T/2, and this is the only case in the 
type of potentials (4.18) that admits all the eigenstates to have a smooth transition to 
the unperturbed ones in the limit c — > 0. One may speculate that the cases where the 
non-zero mode A'-independence is seen, such as the above case where the cancellation of 
the c-dependence occurs, may be found more generally, while those cases where it is not 
seen form a set of measure zero in the total space of the non-zero mode of the potential 
\'{t). 
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5. Conclusion and Discussions 



In this paper we investigated various aspects of classical and quantum caustics for 
quadratic Lagrangians in one dimension, and applied the results thereof to examine Je- 
vicki's correspondence hypothesis (1.3) between stationary points and simple poles in the 
amplitude. 

When caustics occur, classical paths spreading out from an initial point are focused 
into a unique conjugate point; in other words, only those paths connecting the conjugate 
points are allowed classically. The constant fc(A), which is the stretching factor between 
the two conjugate points, is found to be a useful pointer, along with the Morse index 
m(A), for characterizing the intrinsic features of caustics. Correspondingly, in quantum 
mechanics the transition amplitude is nonvanishing only for the conjugate points (for 
which it diverges) . We derived the path-integral kernel for the amplitude in a closed form, 
expressed solely in terms of the stretching factor, the Morse index, and the action of (any 
of the) solution paths. 

Contrary to the situation on caustics, once one goes away from caustics the classical 
and quantum situations become different. Our study focused on how the typical feature 
of caustics, i.e., the concentration of intensity at the conjugate point, can be affected by 
quantum effect. The Gaussian slit (gedanken-)experiment using a Gaussian wave packet 
shows that, although the variance itself is enhanced at the quantum level, the susceptibility 
of the variance of the wave packet against initial momentum fluctuations is suppressed. 
High intensity is realized near caustics of two different types, the first being the caustics 
associated with the period [0, T] while the second being those with [— r, T] which is the 
entire period of our experiment. As a quantum analogue of the Jacobi field, we introduced 
the susceptibility which reduces to the ordinary Jacobi field in the classical limit % — > 
0. The susceptibility gives a measure of stability for the concentration, and our result 
shows that the intensity of the amplitude is stabilized by quantum effect. For practical 
purposes, however, we need to extend our analysis further to the effects of higher order 
terms in the Lagrangian near caustics and thereby supplement earlier works carried out 
on caustics. Another possible direction would be to consider these problems in higher 
dimensions allowing for settings more realistic for physical application. 

The classical and quantum features of caustics studied here have then been employed 
to investigate in detail the role of instantons in quantum mechanics. Viewed as a con- 
strained system on the plane IR 2 , the quantum mechanical rotor over the circle S 1 becomes 
two harmonic oscillators sharing the same, time-dependent harmonic potential A given by 
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the Lagrange multiplier field. Singularities in the amplitude arise at those A under which 
caustics occur. We however learned that generically no simple pole can arise, simply be- 
cause the two stretching factors of the harmonic oscillators do not coincide for generic 
boundary conditions on S 1 , except when the two endpoints are the same or diametrically 
opposite with each other. But even then, the zero-instanton contribution is missing in 
the simple poles to complete the correspondence (1.3). Even if we are confined to the 
former specific cases, and ignoring the missing piece, we come across the problem that it is 
impossible to replace the path-integration over the potentials A by an ordinary integration 
over the constant (zero-mode) part Ao as originally assumed by Jevicki. We observed that, 
in the non-trivial example we looked at in section 4, the assumption of the non-zero mode 
independence of the phase factor does not hold in general, although there are cases where 
it holds thanks to a certain cancellation mechanism. From our result we find it hard to 
sustain the hypothesis even in the toy model, let alone in more general systems such as the 
(LP N model of the original concern. We, however, mention that there remains a possibility 
that, even though it is not strictly true, for some reason the assumption may still lead 
to a correct answer, validating the hypothesis in the end. We feel that this possibility 
deserves a fuller study in order to settle down the issue which is important in exploring 
non-perturbative methods in quantum field theory. 
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Ministry of Education, Science and Culture (No. 09740199). 
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Appendix 



In this Appendix we give a direct argument to show that, for non-critical A, the phase 
(3.19) of the kernel (3.18) is independent of the choice of the function p(t). More explicitly, 
we shall show that the phase is given by the action, 

$(6, a; A) = J cl (6, T; a, 0), (A.l) 

which is evaluated for the classical solution x c \ under a non-critical A obeying the boundary 
condition, x c i(0) = a and x c \(T) = b. Prior to our argument, we recall that the function 
pit) which satisfies p(0) = and p(T) = b — c is introduced to compensate the gap at t = T 
between the given endpoint b and the endpoint c for which an actual classical solution x c \ 
exists. Such a gap is inevitable if A is critical and /c(A) a + s(T) ^ b. However, we stress 
that even for non-critical A it is perfectly legitimate to use the classical solution x c \ with 
x c i(T) = c 7^ b for computing the transition amplitude. Thus, for non-critical A, it must 
be that the apparent p(t)-dependence in the result (3.19) disappear and the total phase 
be given by the classical action for the solution x c \iT) = 6, i.e., (A.l). 

To confirm this, we first consider a solution u{t) of the homogeneous equation (2.3) 
with -u(O) = 0. We may expand it in terms of the orthonormal eigenf unctions {u n } in 
(2.19) as 

u(t) = ^2c n u n (t)+u(T)e(t-T) . (A.2) 

n 

In (A.2) we introduced the step function 6>(t — T), defined by 6*(x) = for x < and 
6>(x) = 1 for x > 0, to account for the condition w(T) ^ 0. Differentiating (A.2) with 
respect to t and setting t = T, we find 

ii(T) =J2 c n <{ T ) + < T ) *(0) • (A.3) 

Note that 



c n = I dtu(t)u n (t) = [ dtu{t) 

JO Jo 



2 



M „ (t) = , (A . 4) 



where we used integration by parts and the equation of motion for u. Combining (A. 4) 
with (A.3) we acquire the formal identity, 

«^ = ff (0) (A.5) 
uiT) ^ E n 
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It is true that the r.h.s. of this identity is ill-defined because of 5(0) and the sum of 
the infinite series, but it may be made sensible by demanding that it be valid for A = 0. 
Indeed, since we have u(t) = at with some constant a for A = 0, the demand suggests 

1 = T 5{Q) - ■ ( A -6) 

n 

This allows us to 'renormalize' (A. 5) to get 

This time the r.h.s. is well-defined, since each of the terms in the infinite series becomes 
small sufficiently fast when n becomes large, as seen from the asymptotic behaviour of the 
eigenfunctions [13, 14]. 

Consider next the expansion of the function p(t) appearing in (3.19) analogous to 
(A.2), 

p(t) = J2 d n Mt) + P(T) 0(t - T) , (A.8) 

n 

where d n = J Q dt p(t)u n (t). Substituting (A.8) in the phase (3.19) we find 

$(&, a; A) = I[x cl ] + x cl (T)p(T) + \p 2 (T) |<J(0) - E ^} • (A- 9 ) 

Now, using the identity (A. 5), and noticing that the solution u may be given by u = 
j3(x c \ — x c \) with (3 a constant, we arrive at 



1 



T 1 



$(6, a; A) = -x cX x cX + -p(T){x cl (T) + x cl (T)} 



2 

1 



02 (A.10) 



Id (b, T; a, 0) + - { a x cl (0) - a x cl (0) + b x cl (T) - c x cX (T) } . 



The relation (A.l) then follows by use of Lemma 1. 
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